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Abstract 

This is an expanded version of f3l. Using techniques from [4 ', in which a differential- 
geometrically intrinsic treatment of counterparts of classical electrostatics was intro- 
duced, it is shown that on some compact manifolds, certain stable configurations of 
points which mutually repel along all interconnecting geodesies become equidistributed 
as the number of points increases. 



In this note, we describe a feature of stable configurations of repelling 
points on certain manifolds in the sense of [4J (pp. 282-288). Our method is 
applicable in several contexts, in particular to flat tori and compact hyperbolic 
manifolds, but we will confine our discussion to the case of compact hyper- 
bolic manifolds, since the detailed relevant background material is already 
available for that case in H. The discussion of this subject in H was mo- 
tivated by analogies with classical electrostatics, but the conclusions drawn 
there can, if preferred, be regarded as being of a formal character, and inde- 
pendent of these analogies. We note that in the flat torus case the analogue 
of the Selberg pre-trace formula is the Poisson summation formula. 

Suppose that AI is a compact hyperbolic manifold, and that is a finite 
set of distinct points in M. We regard the points of S as mutually repelling 
along connecting geodesies, with the repelling magnitude varying with dis- 
tance, as specified by a function H{p) p e [0, oo). In other words, the effect 
exerted on a point p £ Shy points in S is computed by taking the (infinite) 
vector sum H{Ln)Vn, where L„ ranges over the lengths of the geodesic 
segments gn connecting p to points of S, and Vn is the unit tangent vector to 
gn at p, taken in the repelling direction. Note that p is permitted to act on it- 
self. Conceptually, this can be visualized intrinsically on M, or if one prefers, 
on the universal covering space iJ", where the mutual repulsions act along 
the countably numerous geodesies connecting points, attenuated according to 
H{p). Similar considerations can also be applied to continuous, rather than 
point distributions on M, where the repulsion between infinitesimal parts of 
M is defined analogously lIH, but in this note we will be concerned with fi- 
nite point distributions. In the interests of expository economy and to avoid 
unnecessary duphcation, we refer the reader to ||4J for background material 
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concerning the Selberg pretrace and trace formulas, as well as to the ap- 
plication of the pretrace formula to the matters considered here. In greater 
detail, following the discussion in Q, we assume that H{p) is linked to a 
function k{p) by the requirement that H{p) = —k'{p), where k{p), which 
corresponds to a point-pair invariant, is a smooth function on [0, oo) which 
vanishes at infinity, and for which the Selberg pretrace formula is valid, with 
uniform and absolute convergence. As is customary, we denote the Selberg 
transform of k{p) by h{r), and we will be concerned with the values of h{r) 
on the cross-shaped subset of the complex plane composed of the union L of 
the real axis with the closed segment of the imaginary axis from — (n — l)*/^ 
to (n — l)j/2. In the standard parameterization r ^ A, this corresponds to 
eigenvalues of the positive Laplacian on M, which are situated on [0,cx)). 
We will also impose, for the purposes of the present note, the additional re- 
quirement that h{r) is positive on L, and will call the configuration S stable, 
if the net repelling effect at each point of S is zero. 

Let = Ao < Ai < A2 . . . be the eigenvalues of the positive Laplacian on 
M, with repetitions to account for multiplicity, and suppose Lpo,ipi,(p2, ■ ■ ■ 
is an associated orthonormal sequence of eigenfunctions. Denote the points 
of the set 5 by a;i, . . . , xjy. With this notation, it follows from the discussion 
in (|4| (cf. especially p. 287), that sets S for which the quantity 



is locally minimized in the iV-fold Cartesian product of M with itself are 
stable in the above sense. In somewhat more detail, formula (17) from page 
283 of f4\, which describes the repelling effect produced by a point of M 
on another point of M, can be used to derive this result. The requirement 
of zero mutual repulsion at each point of S is expressible via that formula 
as a requirement for the simultaneous vanishing of several instances of (17). 
This translates, since the gradient on the A^-fold Cartesian product of M with 
itself splits into the vector sum of the gradients on the factors, into a necessary 
condition for (1) to have a minimum. 

We remark that there is a misprint on page 283 of |4|. Namely, in 
the formula above (17), the expression —Vxk{x, jy) should be replaced by 
— J2j ^(^; 72/)' ^iid directly above, the words "and by the pretrace for- 
mula" should be augmented to "and by the pretrace formula, counting half 
the r„'s to avoid a factor of i." 

These observations suggest, as being of particular interest, the study of 
the statistical behavior, as — > 00, of point configurations that are glob- 
ally minimizing, and as we will now show, such configurations are equidis- 
tributed in the limit. In order to see this, we will prove a somewhat more 
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general fact, expressed as a theorem, from which the assertion immediately 
follows. 

Theorem. Suppose that M is a compact Riemannian manifold, and that 
l,Lpi,Lp2, ■ ■ - isa sequence of continuous functions inC{M), each having L"^ 

norm 1, and with Jj^ ^Pi = for i = 1,2, Suppose also that finite linear 

combinations of functions from the sequence are sup-norm dense in C{M). 
An important special case in which these conditions are satisfied occurs when 
the sequence 1, (^i, 952, • • • is complete and orthonormal in L^{M). Denote 
by the N-fold Cartesian product of M with itself Let ai, 02, • • • be 
a sequence of positive numbers for which the series Yl^=i Wn{x)\^ an is 
uniformly convergent, and suppose, for N = 1, 2, 3 . . that Xn G 
globally minimizes the quadratic expression 



XI \^n{xi) + . . . + fPniXN)]^ an ■ 



n=l 



Then, as N ^ 00, the sets Xpj become equidistributed, in the sense that for 
any continuous function f on M, 

as N ^ 00. 

Proof. It follows from the hypotheses that for each n, 

/ \^n{Xl) + . . . + ^PniXN)f dXi ■ ■ ■ dXk = iVF^"^ , 

where V is the volume of M, so 
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\^n{xi) + . . . + ipn{.XN)\^ an dxi ■ ■ ■ dxk = NV^ ^ ^ a„ . 
Mx-xM J 

In particular, by the mean value theorem for integrals, there exists {xi, . . . ,xn) & 
such that 

00 TV ^ 

IV^n(^l) + . . . + ^n{XN)\ an= y^^an, 
n— 1 n—1 

which shows that if (pi . . . ,Piv) € is globally minimizing, then 

00 00 

X! \^n{P\) + ■■■ + ^n{PN)\^ an< y On • 
n—1 n—1 

For a fixed m, we therefore conclude from the positivity of the a„'s that 
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so 



\fm{Pl) + • • • + fm{PN)\ < 

which implies that 
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C{m) 



where 



C{m) 



\ 



Va 



^ OO 



In particular, the Xk's, in the limit integrate the uniformly dense set of func- 
tions l,ipi,(p2, ■ ■ ■ (the constant function is automatic), so Weyl's criterion is 
satisfied, and the X^'s are equidistributed in the hmit. □ 
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Remarks. 



1 . It is not difficult to give examples of force laws for which h is positive. 
This is exceptionally easy to see in the hyperbolic 2-dimensional case, 
since then k{u) — ^h{u), so there are many suitable h, k pairs. In 
the flat torus case, k and h are in all dimensions Fourier transforms of 
each other, which again provides many examples. We note also that 
although it is customary to do so, there is no formal reason to insist 
that the force law be always repelling. 

2. If we consider globally minimizing configurations to be equivalent 
when they are mutual images by isometrics of A/, it would be inter- 
esting to know something about the behavior with N of the number of 
inequivalent configurations. 

3. Our definition of point interaction is intrinsic, and does not depend 
on an isometric embedding of M. A prominent example of the latter 
occurs in the well-known Thompson problem concerning stable con- 
figurations of electrons on the unit sphere of R^, where the repelling 
actions resulting from the Coulomb law are transmitted via mutually 
connecting geodesic segments in the containing space R^. Note that 
the Coulomb law has a singularity at the origin, whereas the methods 
of this paper require the function k to have a smooth even extension 
to R^, and the series which arise to be suitably convergent. This of 
course does not rule out the possibility of studying an intrinsic singular 
law as a limit of regular laws. In view of various isometric embedding 
theorems, there can be connections between the intrinsic and extrinsic 
viewpoints. For example, if distance computed in the containing space 
depends smoothly on intrinsic distance, the extrinsic case may in some 
instances become a subcase of the intrinsic case, by choosing the in- 
trinsic repelling function H{p) to smoothly coincide up to the diameter 
of M with the repelling function associated with the embedding, then 
appropriately bringing it down to zero over a distance less than, for ex- 
ample, the injectivity radius of Af, provided this can be done in such 
a way that the positivity requirement on h{r) is satisfied. There can, 
however, be substantial differences between the extrinsic and intrin- 
sic viewpoints. For example, in the case of the unit circle embedded 
in various ways in R'\ stable configurations for the embedded prob- 
lem corresponding to various repelling laws in i?" are sensitive to the 
shape of the embedded circle (cf. [2j), whereas this issue is obviously 
not present for the intrinsic problem. Some asymptotic and extremal 
aspects of the intrinsic problem for geodesically equispaced points on 
a circle are discussed in f5\ and fT\. 

4. In the intrinsic case, it is possible for a point to act non-trivially on 
itself. In more detail, it follows from an examination of the analysis 
leading to formula (17) of [4|, that zero-distance action of a point on 
itself plays no role, but that the effects of the non-trivial geodesies 
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connecting a point to itself are counted. In the flat torus case, these 
actions cancel out in pairs, and so in that case a single point is stable 
with respect to itself, but in the hyperbolic case, this is not generally 
true. In the hyperbolic case, the problem becomes one of studying 
minima for the function 

oo 

5]|^„(x)|^/.(r„) 

n=l 

on M. 
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